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Abstract 

In this paper, an algorithm to compute a certified G 1 rational parametric approximation for algebraic space curves is given by 
extending the local generic position method for solving zero dimensional polynomial equation systems to the case of dimension 
one. By certified, we mean the approximation curve and the original curve have the same topology and their Hausdauff distance is 
smaller than a given precision. Thus, the method also gives a new algorithm to compute the topology for space algebraic curves. 
The main advantage of the algorithm, inhering from the local generic method, is that topology computation and approximation 
for a space curve is directly reduced to the same tasks for two plane curves. In particular, the error bound of the approximation 
^*) ' space curve is obtained from the error bounds of the approximation plane curves explicitly. Nontrivial examples are used to show 
• the effectivity of the method. 
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1. Introduction 

Algebraic space curves have many applications in com- 
puter aided geometric design, computer aided design, and 
geometric modeling. For example, the algebraic space 
curves defined by two quadrics are widely used in geometric 
modeling. One can have an exact parametrization for these 
algebraic space curves. However, exact parametrization 
representation for general algebraic space curves do not 
exist. And usually, we are interested in a rational paramet- 
ric representation. So the use of approximate techniques is 
unavoidable for parametrization of algebraic space curves. 
Some approximate techniques are able to reproduce exact 
rational parameterizations, if those are available. Other- 
wise, one usually approximates an algebraic space curves 
with piecewise rational curves under a given precision. 
Moreover, sometimes, one requires that the approximation 
curves preserve the topology of the original algebraic space 
curve. We call the approximation certified (at precision 
e) if it has the same topology as the original one and the 
Hausdorff distance between the curve and its approxima- 
tion is upper-bounded (by e) simultaneously. 



There are several difficulties for approximate parametriza- 
tion of algebraic space curves. The first one is to preserve 
the topology of the algebraic space curve. In fact, there 
already exist some related work of computing the topol- 
ogy of algebraic space curves, reduced or non-reduced, for 
example [3,16,30,24,33]. Most of them require the curve to 
be in a generic position. For the space curves which are not 
in a generic position, one need to take a coordinate trans- 
formation such that the new space curve is in a generic 
position. Thus some geometric information of the original 
space curve is lost. Some non-singular critical points of the 
new space curve may not correspond to the non-singular 
critical points of the original space curve. One needs ad- 
ditional computation to get these points in the original 
coordinate system. Subdivision method can preserve the 
topology of the curve in a theoretical sense. But it is rather 
difficult to reach the required bound in practice currently 
[11,32]. Even if one gets the topology of the given curve, 
the approximation curve may have different topology as 
the original curve when two or more curve segments are 
very close ( see Figure 2). The second difficulty is the error 
control of the approximation curve. Some error functions is 



reliable but it is not easy to compute in practice, for exam- 
ple [15]. Wc need to find a reliable and efficient method to 
control the error during the approximation. The third one 
is the continuity of the approximation curve. We usually 
require the approximation to be (7 (or G 1 )-continuous or 
higher in practice. Doing so, we need compute the tangent 
directions of algebraic space curve at some points. It is not 
a trivial task especially when the component considered 
is non-reduced. Its tangent direction can not be decided 
by the normal directions (§£, %, |f ), (|f , % ,|f ) of the 
two surfaces / = and g = at the given point. In the 
non-reduced case, two normal directions are parallel or at 
least one of them does not exist at the point. 

There exist nice work about approximating of alge- 
braic space curves. Exact parametrization of the in- 
tersection of algebraic surfaces is obtained in [1,10,17- 
19,22,27,39,41,43,45]. Of course, it is topology preserving. 
The approximation of the intersection of generic alge- 
braic surfaces with numeric method is also considered 
[5,6,31,34,36]. Usually, numeric method cannot guarantee 
the topology of the original algebraic space curve. 

In [7] , the authors considered approximating of the reg- 
ular algebraic space curves with circular arcs by numeric 
method combining with the subdivision method. It works 
well for low degree algebraic space curves. In [23], the au- 
thors present an algorithm to approximate an irreducible 
space curves under a given precision. It based on the fact 
that there exists a birational map between the projection 
curve C for some direction and the irreducible algebraic 
space curve S. But an irreducible decomposition of a given 
two polynomials system is not an easy task. And we need to 
consider the intersection of two or more irreducible space 
curves after we decompose a reducible space curve. Other 
type of intersection problem of surfaces can be found in 
[34,36] and related references. 

In [10], the authors present an algorithm to approxi- 
mate an algebraic space curve, defined by / — g = 0, 
in the generic position with Ferguson's cubic p(t) = 
(x(t), y(t),z(t)),t G [0, 1] and by minimizing an integral to 
control the error. They compute the topology of the space 
curve at first, so it is topology preserving. But they do not 
check whether the approximation curve exactly preserves 
the topology of the original space curve. And it works well 
for regular space curve. From the formula above, we can 
find that if some segments of the algebraic space curve is 
not regular, the method may fail. 

In [37], the authors consider the irreducible algebraic 
space curve in generic position such that its projection is 
birational. They use a genus plane algebraic curve to ap- 
proximate the projection plane curve under a given preci- 
sion if it exists. Thus they have a rational approximation 
space curve for the original space curve. The method is not 
topology preserving. 

In this paper, we present a new algorithm to compute 
a certified G 1 rational parametric approximation for alge- 
braic space curves, which solves the three difficulties men- 
tioned above nicely. The algorithm is certified in the sense 



that the approximation curve and the original curve have 
the same topology and their Hausdauff distance is smaller 
than a given precision. The algorithm works for algebraic 
space curves which need not to be regular or in generic po- 
sitions. The key idea is to extend the local generic position 
method [13,12] for zero-dimensional polynomial systems to 
one-dimensional algebraic space curves. The algorithm con- 
sists of four major steps. 

Firstly, the space curve S, which is the intersection of 
f(x, y, z) = and g(x, y, z) — 0, is projected to the xy- 
plane as a plane curve C\ and C\ is approximated piece- 
wisely with functions of the form x G [a, b\. 

Secondly, we find a number s > such that under the 
shear transformation ip : (x,y,z) — > (x,y + sz,z), </?(/) 
and f(g) are in a generic position in the sense that there 
is a one to one correspondence between the curve segments 
of S and that of their projection curve Ci to the xy-planc. 
The plane curve Ci is also approximated piecewisely with 
functions of the form h^ix), x G [a, b]. 

Thirdly, we choose s such that C2 is in a local generic 
position to C\ in the following sense. 

- The plane curves C\ and C2 can be divided into segments 
such that each segment of C 2 corresponds to a segment 
ofCi. 

- Let h\{x) , fi2(x) , x G [a, b] be the approximations for a 
segment C\ of C\ and the corresponding curve segment 
C2 of C2 with precisions ei and e 2 respectively. Then 
the space curve segment S corresponding to C2 can be 
approximated by (x, h\{x), h ^ x )~ h ^ x ) ^ w ith precision 

^ S 2 £ 2 + (£i + £2)2/s _ 

In other words, if C2 is in a local generic position to C\, 
then each segment of the space curve can be represented 
as a linear combination of corresponding segments of C\ 
and C2. As a consequence, a certified parametrization for 
the space curve can be computed from that of C\ and C2 
directly. This step is the main contribution of the paper. 

Finally, we show that a plane curve can be approximated 
such that the piecewise approximation curve for the space 
curve has G 1 continuity and usually has the following forms: 
(x,ax + b+ J ^ TI ,a 1 x + b 1 + 1 J^ TI + 1 ^ TI ), (ay + b + 
j^,y, a iy + h + + jSL.). 

The topology of the space curve is obtained directly from 
the two projection steps. Thus this new method can not 
only compute the topology of the space curve but approx- 
imate the algebraic space curve under any given precision. 

The paper is organized as below. In the next section, we 
will consider the certified approximation of plane algebraic 
curve under a given precision. In Section 3, we will show 
the theory and algorithm for certified approximation of al- 
gebraic space curves. In Section 4, we will show some ex- 
amples to illustrate the effectivity of our method. We draw 
a conclusion in the last section. 



2. Approximate parametrization of plane algebraic 
curves 

Given a plane algebraic curve defined by a square free 
polynomial / € Q[x, y], our aim is to give a piecewise C 1 - 
continuous approximation of C = {(x,y) € R 2 \f(x,y) = 0} 
in a given box B = {(x,y) £ M. 2 \a < x < b,c < y < d} 
such that each piece of the approximation curve has the 
form (x, h(x)) and the approximation error is bounded by 
a given precision e > 0, where Q, R are the fields of rational 
numbers and real numbers, respectively. And the whole 
approximation curve has the same topology as C. 

2.1. Notations 

In this subsection, we will introduce some notations. 

A point P = (xo,yo) is said to be a singular point on 
C ■ f(x,y) = if f(x ,y Q ) = f x (x Q ,y Q ) = f y (x ,y ) = 0, 
where f x = , f y = ^ . A non-singular point is called a 
regular point. An x-critical (A y-critical) point P = 
(xQ,y ) is a point satisfying f(x ,y ) = f y {x Q ,y Q ) = 
(f(x ,y ) = f x (x ,y ) = 0). So a singular point is both 
a;-critical and y-critical points. The inflexion points or 
flexes of C are its non-singular points satisfying its Hession 
equation H(f) = (see [44]). 

A regular curve segment C of C is a connected part of 
C with two endpoints Po(xo,yo) and P\{xi,y\) (xq ^ x\, 
both Pq, Pi are bounded) and there are no x-critical points, 
y-critical points, flex on C except for Pq, P\. So a regular 
curve segment is convex, monotonous w.r.t. x or y, and 
inside a triangle defined by its endpoints and their tangent 
directions. Let A be the triangle defined by Pq, Pi and their 
tangent lines. An endpoint of a regular curve segment is 
called a vertical tangent point, VT point for short, if 
the regular curve segment has a vertical tangent line at this 
endpoint. 

A parametric curve is said to be C 1 -continuous (G 1 - 
continuous) if the curves are joined and the first deriva- 
tives are continuous (the curves also share a common tan- 
gent direction at the join point). 

2.2. Curve segmentation of a real plane algebraic curve 

In this subsection, we will show how to divide a plane 
curve inside a box B, denoted as Cb, into regular curve seg- 
ments with the form [P (xo,yo), Pi(xi,yi),T (l, k ),Ti(l, kl 
where Pq,P\ are endpoints and Tb,Ti are tangent direc- 
tions at the endpoints. 

We will follow the steps below. 

At first, Compute the topology of Cb- There are 
many related work to solve this problem, such as 
[2,4,8,14,20,28,38,40]. Some methods work well, but they 
need a coordinate system transformation. We prefer the 
methods which do not take a coordinate system transfor- 
mation. 



Second, Compute all the flexes, ir-critical and y-critical 
points of Cb- y-critical points are computed before, x- 
critical points and flexes can be obtained by solving the 
corresponding equations. 

Third, we split the plane curve into regular curve seg- 
ments at cc(y)-critical points or flexes of Cb- An easy way 
to solve the problem is as follows. For all the cc-coordinates 
of x-critical points and flexes, lifting them to split Cb at 
these intersection points. We can find all these endpoints. 

Finally, we represent the tangent direction of any non- 
VT point as (1, fc), k € M\ {+oo, — oo}. It is convenient for 
our approximation. The tangent direction of a VT point is 
defined to be (1, oo). 

Tangent direction computation of singularities. We 

compute the tangent direction of a point close to a singular- 
ity on the regular curve segment to replace the tangent di- 
rection of the singularity. It is easy to compute. Let P(a, /3) 
be a singularity of a planar algebraic curve h(x, y) = and 
C : (x, y(x)),x £ [a, 7] a regular curve segment originating 
from right side (left side is similar) of P. Then the tangent 
direction of C at P is = (1, lim a ._ i . Q + y(x)). In prac- 

tice, we can take some point very close to P on C, which is a 
regular point. Let [a, b] be the isolating interval of a(=^ a, 7^ 
b) . We can use the tangent direction of some regular point 
to replace the tangent direction of C at P. For instance, 
(1, q£' ) = (1, ^) can be regarded as the tangent direction 
of C at P. For the regular curve segments shall the same 
tangent direction, we can take the average value of them as 
their tangent directions, and the center of the isolating box 
of P as the singularity, as shown in Figure 1. When C has 
a vertical tangent direction, t = 00. If | \ > N, for ex- 
ample, N — 100, we can regard the regular curve segment 
have a vertical tangent direction. 

If we cannot distinguish the tangent directions of two 
groups of regular curve segments, we can refine [a, b] to a 
narrower one and recompute the tangent directions again 
until we can distinguish them or they are less than some 
given bounded value r such that \k — k'\ < r, where k, k' 
are tangent directions. 




Fig. 1. Approximate the tangent Fi §- 2 - The approximation curves 
'direction of regular curve seg- chan § e the to P° lo Sy ° f °"g inal 
ments at a singularity ones 

2.3. Approximation of a regular curve segment 

We will give an approximate parametrization of real 
plane algebraic curves. Though Gao and Li [23] have ob- 
tained a rational quadratic approximation of real plane al- 
gebraic curves with B-splines, we need to derive a piecewise 
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approximation curve as (x, h(x)) of real plane algebraic 
curves in order to approximately parameterize real space 
algebraic curves in a different way. Let C be the regular 
curve segment defined by two points Pq (xo > 2/0 ) > -Pi (xi , yi) ■ 
We divide the approximation problem into two cases. 
C not containing a VT point. Given a regular curve seg- 
ment C, its two endpoints Po(xo, Uo), Pl(xi, TJi) does not 
contain a VT point. And the tangent directions of the regu- 
lar curve segment at Pq, Pi are (1, feo), (1, fci), respectively. 
We will construct an explicit rational quadratic function 
Y\(x) to approximate C such that Yi(x,) — yi,Y{(xi) 
ki, i = 0,1, where Y{(x) 
Assuming that 



Yi(x) 



ax 2 + bx + c 
dx + 1 



(1) 



and xo — 0, x\ — 1, we have 

Ki(0) = y ,^i(l) = yiX(0) = koX(l) = h. 

Solving a, b, c, d from the equations above, we have 

a= -2yoVi +yo 2 - kpkj + yi 2 

-yi +ki +y 
, -2 y y t + y k t + 2 y 2 + k y t - k k t 



c = yo, 
d = -~- 



-yi + kt + y 
-2 yj + k + kj + 2 y 



-yi +kt + y 

From the representation, we need to require 

-Vi + k 1 +y o ^0, (2) 

and dx + 1 has no roots in [0, 1], that is, d > —1, from 
which we can derive that 



(-y t + k + y ) {-yi + k t + y ) < 0. 



(3) 



From the mean value theorem, we know that y\ — yo = 
k x , where k x is the tangent direction of some x G (0, 1). 
Since C is monotonous, so k x is some value between kg and 
fci. Thus conditions (2) and (3) are satisfied directly. We 
can easily transform the interval [0, 1] to [x , xi] by setting 



X — TO 



x = Z" , where x <G [0, 1] when X € [xq,xi]. 



Furthermore, when d ^ 0, that is to say —2y 1 +k + 
kj + 2 y 7^ 0. Then expression (f ) can be transformed into 



Yi (a;) = a x + b + ^ 



dx + 1 



(4) 



Though equation (4) is equivalent to equation (f) when 
d ^ essentially, it has a simpler form and can reduce 
computation when evaluation. When d — 0, equation (1) is 
a polynomial of degree two. And we have simple expressions 
for parameters a, b, c, that is a — k '~ ko , b = k , c = yo. 
C containing a VT point. When a given regular curve 
segment contains a VT point, it means the tangent line at 
Pq or Pi is a vertical line x — x = or x — x\ = 0. In 
this case, the method above does not work. But we can use 



part of an ellipse or a hyperbola 



(x-x a f 



± 



(V-Vaf 

b 2 



- 1 = 



0(a > 0, b > 0) to derive an approximate parametrization 
of a real plane algebraic curve. Note that a regular curve 
segment containing a VT point has four cases which exactly 
correspond to the four parts of an ellipse or a hyperbola: 
the vertical line is x — xq = or x — x\ =0 and y > 
y or y < y a . We consider the case that C has a vertical 
tangent line at Po(^o,2/o) and C monotonously increases 
from Pq to Pi . And we assume that the tangent line at Pi 
is ki(x — xi) — (y — j/i) = 0(fci > 0). The approximate 

curve is Y 2 {x) = y a + b "^ a p ■ . Note that we have 
x a = xq ± a, j/o = yo from the property of the ellipse or the 
hyperbola. So we have 



Y 2 {x) =y + 



b \/\a 2 - (x - x - a) 2 \ 



(5) 



And Y 2 (xi) = yi,y" 2 '(xi) = k\. Solving it, we have 

(xi ~ xo) (xo ki + yi - yo - xi ki) 

a = , 

yi - yo -2xi ki +2 xo ki 

b = (xo ki +yi — yo - xi ki) 



yi - yo 



yi - yo —2xi ki +2 xo ki 



From the representation, we can find that a, b are well 
defined if ki < 2 fx~-x ) ^ or an cm P se or &i > 2(x7-x Q ) 
for a hyperbola. So we can choose {x\,y\) on the regular 

Then we can use 



curve segment such that ki + „ f 1 Vo , . 
part of an ellipse or a hyperbola to approximate the regular 
curve segments with VT points. The other three cases can 
be solved in a similar way. 

Lemma 1 The two kinds of approximation curves above, 
say C : (x, h{x)),x <G [xo,xi] for C are inside the triangle 
formed by the endpoints and the tangent directions at the 
endpoints of C , denoted as A. 

Proof. For both cases, C is part of a quadric curve. And the 
curve intersects all three edges of A at least twice (including 
the multiplicities). If C goes out of A, it will intersect the 
edge(s) at least three times (including multiplicities). But 
it is not possible since C is part of a quadric curve. So the 
lemma is true. □ 

Topology preserving approximation. After we get the 
approximation regular curve segments, we need to check 
whether the approximation curve change the topology of 
the original curve. Even if we get the correct topology of 
the given algebraic planar curve, the approximation curve 
may have a different topology as the original curve, espe- 
cially when two regular curve segments are very close, for 
example, see Figure 2. So we need to ensure that our nu- 
meric approximation curve has the same topology as the 
original one. We need only ensure that any two approxi- 
mation curves, say C\,C 2 : (x,p(x)),(x,q(x)),x <E [a, b], 
are disjoint. If p(x) — q(x) = has no real roots in (a, b), 
then the two approximation regular curve segments are dis- 
joint. There are two kinds of approximation curves, say 
Yi(x), Y 2 {x). So we need to consider: 

Case one: two approximation curves are both rational ones 
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as Y\(x). Then T(x) = p{x) — q(x) can be simplified into a 
cubic univariate polynomial. It is easy to check whether it 
contains a real roots in (a, b) by its coefficients. 
Case two: one is as Yi(x) and the other is as Y 2 {x). Then 
T(x) = p(x) — q(x) can be simplified into a quartic univari- 
ate polynomial. It is also easy to check whether it contains 
a real roots in (a, b) by its coefficients. 
Case three: both approximation curves are as Y 2 (x) . They 
both arc parts of quadric algebraic curves. Considering the 
intersection of two quadric algebraic curves, we can judge 
whether Ci, C 2 are disjoint or not. 

Doing so as above, our approximation is exactly topology 
preserving. 

2.4. Error control of the approximation 

We will show the error control of the plane approximation 
curve in this subsection. In geometry, the approximation 
error should be defined as the following Hausdorff distance 
between the segment S and its approximation S a , 

e(S, S a ) = dis(S, S a ) = max^min d(P, P'). (6) 

However such a distance is difficult to compute. As an im- 
plement, the distance from an approximation parametric 
curve P(t) = (x(t),y(t)),0 < t < 1 to the implicit defined 
curve C : f(x, y) = is taken in the following form, which 
is called the error function [15], 

e(t) = 7(x(t),^)) (7) 

^{f x (x(t),y{t)? + f y {x{t) lV (t)) 2 ' 

The approximation error between P{t) and C is set as an 
optimization problem 

e(P(t),C) = max (e(t)). 

0<t<l 

Let C : (x,y(x)),x e [xo,Xi] be the regular curve seg- 
ment and C := (x,Y(x)),x G [xo,xi] its approximation 
curve. It is not difficult to find that the following bound 
is an upper bound of the Hausdorff distance between the 
segment C and its approximation curve C from (6): 

max \Y(x)-y(x)\. (8) 

xE[xo,xi\ 

We use Newton-Ralphson method to obtain y(x®) at 
some point x° € [xo,Xi] in practice and Y{x { () is the start 
point. If we fail to get a point with Newton-Ralphson 
method or the point satisfying |F(x°) — y(xf)\ > 5, we 
can divide the regular curve segment into two ones. The 
approximation error is bounded by maxi {\Y(xf) — y(xf)\} . 
In practice, we sample x° as x? = x + i/n (xi — x ), < 
i < n, for a proper value of n. 

In order to control the error under a given precision, we 
need to divide the regular curve segment into two or more 
regular curve segments recursively until the error require- 
ment satisfied. We subdivide the regular curve segments 
into two or more regular curve segments uniformly in the x 



coordinate. For any regular curve segment C : (x, y(x)),x G 
[a, b], we denote the endpoints as Po(xo, yo), P\{x\, y\) and 
the tangent directions as (1, i = 0,1. We can find that 
Po, Pi and two tangent directions form a triangle. One can 
subdivide the regular curve segment into two (or more) 
ones, for example, C\ : x € [xo,(xo + Xi)/2],C2 : x € 
[(xo + xi)/2,xi], if the precision is not satisfied. For the 
approximation curve of our method, we will prove that it 
can achieve any given precision. 

Theorem 2 Let Pq , Pi be the endpoints of a regular curve 
segment C and A the triangle related to the regular curve 
segment as defined before. The Hausdorff distance between C 
and its approximation curve(s) tends to zero if we subdivide 
C into two or more regular curve segments recursively. 
Proof. We will consider divide C into two regular curve 
segments for the proof since dividing them into more reg- 
ular curve segments are similar. Let P((xo + xi)/2,y) be 
a point on C. Denote the triangles formed by Pi, P(P, P2) 
and the tangent directions of C at these points as Ai(A2). 
Let the lengths of the line segments PiP,PPi be L\,L\ 
and the heights of the triangles Ai, A2 corresponding to 
the edges PiP, PP2 are H\,H\, as shown in Figure 3. Sub- 
dividing the regular curve segments recursively in a similar 
way, denoting the length of the edges and heights as Lj , Hj 
of the triangles, we have the sum of the areas of these tri- 
angles are 

3 3 

Assume that one edge L l - = PfP^, Pi(x' 1; y[), P^x'2, y' 2 )- 
Since the given regular curve segment is bounded, \y' 2 — y[\ 
tends to zero when \x' 2 — x\\ tends to zero. And ^ . L % - 
tends to the length of arc, say £, of the given regular curve 
segment and H % - tends to zero when all corresponding \x' 2 — 
x'i\ of Lj tends to zero. Thus A tends to zero. From the 
result of Lemma 1 , we can have the opinion that we have a 
proof of the theorem. □ 

There are two ways to find the subdivision points on a 
given regular curve segment C : (x, y(x)),x € [a, b]. Since 
we get the topology of the plane projection curve, we know 
the order of the given regular curve segment among all the 
regular curve segments of the projection curve h(x, y) = 
when x changes from a to b. That is, we can find the point 
on C for a fixed x coordinate, say xo € (a, b). It is the real 
root with the same order of h(xo,y) — in a fixed interval 
(or (-00, +00)). 

Another way is a local method. We can trace the regular 
curve segment to find the point on C with given x coordi- 
nate since the regular curve segments are monotonous and 
convex. From the endpoint of the regular curve segment, 
compute the tangent line of the regular curve segment at 
some point P, find a point Q on the tangent line by increas- 
ing the x coordinate such that the line segment PQ has no 
intersection with the projection curve. Then fix the x coor- 
dinate of Q, to find a point on the regular curve segment. 
Note that we know the direction to find the point from the 



positiveness or negativeness of the tangent direction. Do- 
ing so recursively, we can find the point that we want, as 
shown in Figure 4. 




Fig. 3. Splitting a regular curve Fig. 4. Finding subdivision point 
segment by tracing 

With the preparation above, we have the following algo- 
rithm to approximate a plane algebraic curve. 
Algorithm 3 The inputs are C : f(x,y) = 0, a bounding 
box B and an error bound 5 > 0. The outputs are para- 
metric curves C \ := {Bi(x) = (x,yi(x)),a.i < x <bi,(i = 
1 , . .., N) }, such that they give aC 1 - continuous and topology 
preserving approximation toCs with e(C\,C) < S. 

(i) Regular curve segmentation of Cb as in Section 2.2. 

(ii) Regular curve segment approximation as in Section 
2.3 with error control of the approximation as in Sec- 
tion 2.4. 

The correctness of the algorithm is clear from the analysis 
above. The termination of the algorithm is guaranteed by 
Lemma 1 and Theorem 2. 

3. Certified approximate parametrization of 
algebraic space curves 

In this section, we will consider approximate parametriza- 
tion of algebraic space curves defined by /, g £ Q[x,j/, z] 
such that two assumptions hold: 

- For any xq £ K, f(xo, y, z) = g(xo, y, z) = has a finite 
number of solutions; and 

- the leading coefficients of /, g w.r.t. z have no common 
factors only in x. 

The assumptions are to ensure that we can use local generic 
position method to recover the points on the space curve 
from the points on two plane projection curves. The first 
assumption ensures that the algebraic space curve defined 
by / = 9 = does not have a plane curve on the plane 
x — xq. The second assumption ensures that we can find 
a generic position only by taking a shear map (x, y, z) — > 
(x, y + s z, z). If the two projection curves have no factors 
only involving x, the two assumptions hold. 

In fact, most of the problems we considered satisfy the 
condition. Note that we can exchange x, y, z freely. And 
another coordinate system transformation (x, y, z) — > (x + 
s z, y, z) can help us to find out the missing regular curve 
segments in the first transformation, even when the alge- 
braic space curve containing vertical lines. Thus we can 
remove the assumptions with the method mentioned here. 
But we still assume the two assumptions holds in this sec- 
tion. 



3.1. Definition of local generic position 

In order to reduce the 3D approximation of space curves 
into 2D approximation of plane curves, we need the concept 
of local generic position. We recall the related definitions 
for zero-dimensional bivariate polynomial system [12]. Let 
C be the field of complex numbers. Let /, g £ Q[x,y]. We 




Fig. 5. LGP method 

say two plane curves defined by two polynomials /, g such 
that gcd(/, g) — 1 are in a generic position w.r.t. y if 

1) The leading coefficients of / and g w.r.t. y have no 
common factors. 

2) Let h be the resultant of / and g w.r.t. y. For any 
a £ C such that h(a) = 0, f(a,y),g(a,y) have only one 
common zero in C. 

Then we will introduce the technique of local generic 
position (LGP for short) method. 

Given f,g£ Q[x, y], not necessarily to be in generic po- 
sition, we can take a coordinate system transformation 4> : 
(x,y) —> (x + sy,y),s £ Q such that 

- </>(/), 4>(g) are in a generic position w.r.t. x. 

- Let h, h be the resultants of <fi(f),4>(g) and f,g w.r.t. 
y, respectively. Each root a of h(x) = has a neighbor 
interval H a such that H a n Hp — for roots (3 ^ a of 
h = 0. And any root (7,77) of / = g = which has a 
same x-coordinate 7, is mapped to 7' = 7 + sn £ H~ n 
where ^1(7) = 0, ^(7') = 0, as shown in Figure 5. Thus 
we can recover n — 7 ~ 7 ■ 

We can find the method has two nice properties: 1) The 2D 
solving problem is transformed into a ID solving problem. 
2) The error control of the solutions is easier. 

3.2. Basic idea 

Now we want to extend this technique to 3-D case. Let 
/ A g denote the algebraic space curve defined by / = g = 
0. Denote ir z : (x,y,z) — > (x,y) and h = Tr z (f A g). Let 
<P ■ (x, y, z) -t (x,y + s z, z), h = Tr z ((p(f) A <p(g)). For a 
proper regular curve segment C : (x,p(x)),x £ [xq,Xx] of 
the plane curve defined by h — 0, it corresponds to one 
(the corresponding space regular curve segment may be at 
infinity) or more space regular curve segment (s), denoted as 
Si, . . . , St- If we can choose a proper s such that <£>(/), <f(g) 



is in "a good" position, and it has some local property, that 
is, the corresponding projection regular curve segments of 
Si,..., S t , say Ci : {x,q t (x)),x <E [x ,xi](i = l,...,t), 
are in a fixed neighborhood of C, then we can recover the 
z-coordinate of St : z = 2il£iz£M(j = ].,...,£). But to 
make C'i in a fixed neighborhood of C is an unreachable 
task sometimes. Fortunately, what we need is to find out 
the correspondence between the regular curve segments of 
the plane curves defined by h = and h = 0. We can 
choose some sample points on x-axis, say x^ ,i — 1, . . . , n, 
a proper s such that (p{f(x^\y, z)),ip(g(xf ) \y, z)) arc in a 
generic position. Then we can figure out the correspondence 
between the regular curve segments of h = and h = 0, as 
shown in Figure 6. 

To realize the aim above, there are two key steps. One 
is how to find an approximation s, the other is how to find 
the correspondence between the regular curve segments of 
h = and h = 0. We need some preparations at first. 

We say that two algebraic surfaces defined by /, g e 
Q[x, y, z] such that gcd(/, g) = 1 are in a z-generic posi- 
tion if 

1) The leading coefficients of / and g w.r.t. z have no 
common factors. 

2) Let h be the resultant of / and g w.r.t. z. There are 
only a finite number of zeros (a, /3) G C 2 such that (a, (3) is 
not a y-critical point, h(a, ft) = and f(ce, ft, z), g(a, f3, z) 
have more than one distinct common zeros in C. 

The definition is similar to the definition of pseudo- 
generic position in [16]. In Theorem 4 of [16], the authors 
also provide a method to check whether two given surfaces 
are in a pseudo-generic position or not. 
Computing s. We will show how to find an s mentioned be- 
fore. Let 7r y : (x, y) — > (x). Denote the real roots of ff v {h) = 
and the x-coordinates of the flexes and x-critical points 
of h = as oil, ots, ■ ■ ■ , a 2 t-i- Find two rational numbers 
less than a\ and larger than a2t-i, denoted as a>o,a 2 t re- 
spectively. For any two adjacent real roots a 2i -i,a 2i+ i of 
n y (h) = 0, we can find a rational number, say a 2 i. Then we 
obtain a sequence ai(i = 0, . . . , 2 1). Assume the real roots 
of h(ai,y) = arc f3i,j(j — 0, . . . ,ti) which are listed in 
increasing order. We can find out that (ai,f3i.j) divide the 
plane curve h = in the region [a , a 2t ] x R into regular 
curve segments. Let 

fi= max RB,(/( aj ,ft,^)), 

0<i<2t,0<j<t i 

r = min {R, min (ft,j + l - ft,,')}, 
0<K2( 0<j<t i -l 

< s < —,s S Q, (9) 
2 R 

where Pi-i = — oo, RB z (/(ai, (3 i:J ;, z)) is the root bound 
of f{oti, fii.j, z) in z, f can be replaced by g. Since it 
is probability 1 under condition (9) to obtain such an 
s that f(f(ai,y,z)),ip(g(ai,y,z)) are in a generic po- 
sition with the assumptions, it is probability 1 that 
<p(f(ai,y,z)),(p(g(ai,y,z)) are in a local generic position 
for all on (i — 0, . . . , 2 1). And we can ensure this by check- 
ing whether <p(f) A ip{g) is in a z-generic position. 



Finding the correspondence. With local generic posi- 
tion method and the assumptions, we can recover the points 
of / A g corresponding to (c^, fiij), say (a i5 fi it j, 7i,j,fe)(l < 
k < tij). As shown in Figure 6, from B(a, (3), Bi(a, (3i), 
we can find out the point corresponding to B,B\ in 3D 
space: (a, P, Note that B l is in a neighborhood 

a x (fi-rjlj- r/2) of B. 

Let a, a' be any ai,ai+i. We will classify the piece of 
curves inside (a, a') x R into two cases by considering 
whether they contain singularities. 

If the two endpoints of a regular curve segment C of h = 
are in the fixed neighborhood of the endpoints of a regular 
curve segment C of h = respectively, we know that C 
corresponds to C, see RiUi(i = 1,2) and RU in Figure 6 
for example. 

There are two cases for the singular points of h = 
w z ((p(f) A tp(g)) — in (a, a') x R: One case is that some 
correspond to singularities of / A g. Two or more space 
regular curve segments of / A g intersect on a cylinder 
surface defined by some factor (s) of h = 0. So this kind of 
singularities of / A g may not correspond to singularities of 
h — 0. If two or more left (right) branches of a singularity 
of h = correspond to a same regular curve segment of 
h = 0, we can judge that it is a true singularity of / A g, 
see the point D in Figure 6 for example. The regular curve 
segments AiDL 2 ^A 2 DLi, A3EFL3 belonging to h — all 
correspond to AR since Ai, A 2 , A^^Ri, R 2 , R3*) are in a 
neighborhood of A{R). 

The other case is that they are not true singularities of 
fAg. Thus the curve branches of h = which pass through 
these singular points of h — with different tangent lines 
correspond to disjoint space curves of <p(f) A <p(g), then we 
can find out the correspondence of C and Ci , see points E, F 
in Figure 6 for example. Note that the continuous space 
curve maps to a continuous plane curve by 7r z . If there exist 
two or more curve branches having same tangent lines at a 
singular point of h = in (a, a') x R, see point G in Fig- 
ure 6 for example, we call it tangent false singularity. 
If Ci contains only one tangent false singularity in (a, a'), 
we can still find out the correspondence of Cj and its corre- 
sponding regular curve segment C following the correspon- 
dence of the endpoints of C and d. Note that if the end- 
points of two regular curve segments of fAg have same x, y 
coordinates, their corresponding projection curves overlap 
in h — 0, and the projection curve of their corresponding 
regular curve segment in </?(/) A ip(g) are disjoint in h = 
(except the endpoints). For G in Figure 6, we know the 

endpoints of_R 3 GC/3 are in the fixedjreighborhood of R 
and U, so R3G, GU3 correspond to RU. But if Ci contains 
two or more tangent false singularities in (a, a'), we can 
not determine the correspondence of the part(s) of Ci be- 
tween these tangent false singularities and C (or other reg- 
ular curve segment of h — 0). As shown in Figure 6, H, K 
are two tangent false singularities and we do not know the 
correspondence of the two regular curve segments between 
them. We can find back the correspondence in the following 



1 



way. Let (po,Pi) and (go, Qi) be two adjacent tangent false 
singularities on Ci in (a, a'). Choose a rational number 7 
such that po < 7 < q Q . Solving f(-f, y, z) = 5(7, y, z) = 0, 
we can get some real points on / A g. Solving h(j, y) = 0, 
we can get some real points on h — 0. Since ip(f) A <p(g) 
is in a z-generic position, ^(/(7, y, z)) A ip(g(j, y, z)) is in 
a generic position. So two group of points have a one-to- 
one map. Thus we can find out the correspondence between 
them. So we can decide the correspondence between C and 
(parts of) Ci. 




Fig. 6. Two projection curves and their correspondences 

Based on the analysis above, we can write a theorem be- 
low which tell us an algorithm to obtain the correspondence 
between the two projection curves h = and h = 0. 
Theorem 4 Use the same notations as above. There exists 
an algorithm to find a proper s to obtain h = 0, and to find 
out the correspondence between the regular curve segments 
of h = and h = . 
Some remarks for the theorem: 

(i) If we consider the algebraic space curve inside a box 
B = [ai, 61] x [02, fa] x [03, fa], we can consider only 
the points, including the x(y)-critical points and 
flexes inside B. We also need to decide the intersec- 
tion between the space curve and the boundaries of 
the box. We can use two planes z — a 3 = 0, z— fa = 
to intersect /, g and project the intersection into xy- 
plane, respectively. Replace h by the product of these 
projections and h. Use two lines y — a-i = Q,y — fa = 
to intersect h = 0, we get some of the projections of 
the boundaries of the algebraic space curve inside B. 

(ii) In fact, this theorem gives an algorithm to compute 
the topology of an algebraic space curve. 

(iii) In the computation in practice, (oti,/3ij) are rep- 
resented by isolating intervals. The corresponding 
method to get r, R can be found in [12]. 

If we approximate the two plane projection curves with 
the forms (1), (5), we can get the piecewise approximation 
parametric space curves of / A g from the correspondence 
of the regular curve segments in the two plane curves in 
Theorem 4. They have the form: 



, , v 2/2(2;) -yi(x) 

(x,yi{x), ),x € [a,b\ 

s 

where yi(x), y%(x) are as Forms (1), (5). 



(10) 



We need also check whether our approximation space 
curve changes the topology of original space curve or not. 
Since the plane approximation curve does not change the 
topology of the plane projection curve, we need only to 
check whether two approximation space regular curve seg- 
ments having the same y coordinate are disjoint or not. 
We assume that the two approximation space regular curve 
segments are 

„ , , s V2(x)-y(x) yi(x)-y(x) 

Ci : (x,y(x), ),C 2 : (x,y{xj, ),x e [a,b]. 



They are disjoint if V2 ^ v ^ — v ^ = has no real 

roots in (a, b). We use the similar method as we check two 
plane approximation regular curve segments in Section 2.3. 



3.3. Error control of the approximation space curves 

In this subsection, we will consider how to control the 
error of the approximation space curve. 
Theorem 5 Use the notations as before. If we approxi- 
mate the plane curves h = and h = with errors ei, ea> 
respectively, the error of each coordinate of the approximat- 
ing curve of the algebraic space curve f A g is bounded by 
max(e 1; ei ^ £2 ), and the Housdorff distance error of the ap- 
proximating curve is bounded by ^ s e i + ^ ei+£2 ) _ 
Proof. Let C : (x,y(x)) (Ci : (x,y~i(x))) be the regular 
curve segment of h = (h = 0) and C : (x,p(x)) (Ci : 
(x,q(x))) its approximation curve, x <G [xo,a;i]. Ci corre- 
sponds to C. Let Si : (x, y(x), z(x)) (exact representation) 
and Si : (x,p(x), ^ x g ^ Xq ^ a S p ace regular 

curve segment and its approximation. From the condition, 
we have e(C, C) < ei,e(Ci, Ci) < e 2 . The error here is de- 
fined by (8). Let us consider the three coordinates of one 
part of the approximation curve S : (x,p(x), q{x) s p{x) ),x € 
[x ,xi]. 

The errors of the first and second coordinates are and 
ei, respectively. For the third coordinate, we have z(x) = 

iii^tzm. Thus 



,., , «(*) - pW , fi(«)-v(») j(>)-p(») 
kw 1 = I 

s s s 

< U/iO) - gWj + \y( x ) - pO)I < ei + £ 2 



So the third coordinate is bounded by ei+ s t2 from (8). 
From the definition of Hausdorff distance (6), we have the 
Hausdorff distance of Si and Sj: 



e(Si,Si)= max min d(P,p')< max d(P, P 1 ') 



y/q + + e2 )A) 2 



v / s 2 e 2 + ( ei+e2 )2 



This ends the proof. □ 

If the required precision for the approximation curve is e, 
we can approximate the plane algebraic curves h = and 



h = with precision 



e from the theorem. 



3.4. G 1 -continuous rational approximation space curve 

We will derive approximation space curve from plane 
approximation curve. And we will re-parameterize the 
non-rational parametric curve into rational ones. Thus 
the obtained approximation space parametric curves are 
G 1 -continuous and rational. 

Lemma 6 Use the notations as before. If we approximate 
the plane curves h = and h = with C 1 -continuous 
parametric curve, the approximation curve of the algebraic 
space curve f A g is G 1 -continuous. 

Proof. Let (x,p(x)),(x,q(x)),x € [x ,xi] be two cor- 
responding approximation curves of the regular curve 
segments of ft = and ft = and p{x),q{x) are 
C 1 -continuous in [xo,xi]. We can obtain the approxi- 
mating curve of the space regular curve segment: S : 
(x,p(x), iMz£M ) ; x € [xo,Xi]. The tangent direction of 
S at any x is (1, (|| — || )/s). From the definition of 
C 1 -continuous, we can find that S is C 1 -continuous since 
(x,p(x)),(x,q(x)),x € [xo,xi] is C 1 -continuous. For the 
3D point P of / A g corresponding to a VT point, if we 
require the approximating space curve is C 1 -continuous 
at P, then the whole approximating space curve is also 
C 1 -continuous. □ 

When re-parameterizing the approximation space reg- 
ular curve segments into rational ones, we need to know 
the tangent directions of the endpoints of space regular 
curve segments. For the endpoints corresponding to non- 
VT points, we can directly get it from the tangent direc- 
tions of the plane curves. For the endpoints corresponding 
to VT points, we can get the tangent directions as follows. 
At first, we assume that (x,p(x)), (x, g(x)),x € [x ,xi] are 
parametric plane regular curve segments of exact algebraic 
regular curve segments (x,y\(x)),(x,y2(x)),x e [xo,xi] 
and xo corresponds to a VT point. The exact tangent di- 
rection of the algebraic space regular curve segment at xo 
is (1, y[(xo), ^ 2 ^ x °^ zMii^°l from the parametric representa- 
tion. Note that (1, oo) corresponds to (0, 1) for plane regu- 
lar curve segments. So for the approximation tangent direc- 



tionatx :(l,^,(^ 



^o))/ S ), if ^2 is larger 



dp(xp) / dq(x ) 
dx ' v dx 

than (or less than) some given value, for example, 100 (or 
-100), we can reset the tangent direction as (0, 1, ( dq g^ — 

^ £ t 1 )/(s ^t 1 ))- Moreover, if (^1 _ ^l)/(s %*) 
is larger than (or less than) some given value, we can set the 
tangent direction as (0, 0, ±1). So the tangent directions at 
x is as (0, l,p),p 7^ or (0, 0, ±1). 

Reparametrization of space curve. If the tangent di- 
rection at xo is (0, l,p), we can re-parameterize the space 
curve segment with the form 



_ , ni f + bi t + ci a 2 t 2 +b 2 t- 

P(t) — ( ; ■ — , t, - 



di t - 



d 2 t + l 



d 3 t- 



-), t e [o,i], 



(ii) 



such that it is G 1 -continuous with other regular curve 
segments at the endpoints. Assume that the two endpoints 
are (xi,yi, Zi), i = 0, 1 and the given tangent directions at 
two endpoints are (x[, y[, z[), i = 0, 1. Thus x' = 0. Here 



for simplicity, we assume that yo = 0, y\ = 1 since we can 
set t = y ^ yo . Bisecting the regular curve segment ensures 



that y[ since the regular curve segment is monotonous. 

We require that the parametric space curve satisfying G° 
and G 1 conditions at the two endpoints. So we have eight 
valid equations from the following equations. 



p Mlt=o = (zo, vo, zo), p (*)lt=i = Oi> vi, 

9P(t) SP(t) 1 , , , 

— — — \t=o = (0, 1, p), — — — |t = i = —t(.*i,Vi, 
dt dt y. 

Solving them, we have one solution as below. 

xq 2 — 2x( f xi + x\ 2 



-xi + X + XQ 



1 / / I r I 

a 2 = (zi<i3 + z 1 d 2 d 3 - z ± d 3 + d s z + z + z 1 d 2 - z 1 d 2 + d 2 zq + z 1 

-2z! + 2z ), 

xq ^2xq — 2 xx + x'^j 



bi 
b 2 



ci 
C2 



— xx + x' x + xq 

1 / / 2 

: -( — z\d 2 + d 2 ZQ + z x d 2 — 2 zxd 2 d 3 + 2 d 2 ZQd 3 + z x d 2 d 3 

d 3 

+2z[d 2 d 3 — zxd 2 d 3 2 + 4 d 3 z Q + 2 zq — 2 z\ + z' + z[d 3 2 + 2 z[d 3 
-4 zid 3 + 2 d 3 z' + z[ - 2 zxd 3 2 + 2 d 3 2 z ), 

■ xo, 

1 i 1 i i i 

rr, r( -2 z l d 2d3 + z 1 d 2 d 3 + 2 z x d 2 d 3 + z x + z a - 2 zi 

d 3 2 (-d 3 + d 2 ) 

+2 zq - z x d 2 - 4 z ± d 3 + z[d 2 + z[d 3 2 + 2 z[d 3 - zxd 2 d 3 2 + 2 d 2 z d 3 

+d 2 z + 4 d 3 z + 2 d 3 z + 2d 3 z + d 3 z - 2 z 1 d 3 + d 3 zq), 

1 ,2111 

■- — — "( _2 z l d 2d 3 + z 1 d 2 d 3 + 2 z 1 d 2 d 3 + z 1 + z - 2 z ± 

d 3 2 (-d 3 + d 2 ) 

+ 2 zq — z\d 2 — 4 21CZ3 + z[d 2 + z' x d 3 2 + 2 2^3 — 2!d2^3 2 + ^2^3 2 20 

+2 d 2 zo< 1 3 + d 2 2 + 4 d 3 z + 2 d 3 2 2 + 2 d 3 2 + d 3 2 2 - 2 Z!d 3 2 ), 
2 xq — 2 xi + x{ 



-xx + x + xq 



where di , are free variables. At first, we require that x\ — 
%i + xo 7^ 0, ds(d,2 — ds) 7^ since they are denominators. 
Second, we require that di t + 1 = 0(i = 1, 2, 3) in t has no 
root in [0,1], that is, di > —1. For i = 1, we have equal 
conditions: 



(x - xi)(x - xi + Xi) < 0. 



(12) 



Since the given planar regular curve segment is monotonous 
(w.r.t. both x and y), the first condition and (12) hold 
directly. We can choose proper di, d^ such that conditions 
hold. 

For the case of tangent direction is (0,0, 1), we can set 
the parametric regular curve segment as 



p(t) = ( 



ai t + 6i t + ci a 2 t + b 2 t + ( 



d 2 i + 1 



d 3 t + 1 



,*),*£ [0,1], 



and solve a similar equation system to get the parametric 
regular curve segments. 

The left problem is to control the precision. Let e be 
the required precision for the whole approximation para- 
metric curve. If the non-rational parametric curve <Si : 
(x,p(x),q(x)),x € [xo,xi] approximates the regular curve 
segment of algebraic space curve S with precision e/2, and 



n 



the new rational parametric curve S2 ■ (x(t), y(t), z(t)) ap- 
proximate Si : {x,p(x),q{x)) with precision e/2, then 1S2 
approximate S with precision e. 

We need to control the approximation precision of S2 to 
cSi . In [42] , the authors consider the approximation of 3-D 
parametric curve with rational Bezier curves. For our prob- 
lem, we need rational curve. For any fixed x° G [xq, x{\, we 
can derive a univariate polynomial equation in t of degree 2 
by p(x ) = y(t). Solving it, we have two real solutions (the 
solutions do exist). Choose the one such that x(t) close to 
x°, say t°. Denote the distance between (x°,p(x°),q(x )) 
and (x(t°),y(t°), z(t )) as D(xq). From the dehnition (6), 
we can find that ma,x x o e ^ Xg Xl ^ D(x°) > e(iSi,52) is an up- 
per bound of the Hausdorff distance of 1S1 and 1S2. We can 
choose some sample points to estimate the error between 
1S1 and 1S2 ■ 

Thus, in the end, we get a G^continuous piecewise ra- 
tional approximation space curve under a given precision. 

When we approximate a regular curve segment contain- 
ing a VT point in practice, we usually select a short dis- 
tance for it since the error control is much easier. 



4. Algorithm and examples 

In this section, we will give the main algorithm to ap- 
proximate algebraic space curves and use some non-trivial 
examples to illustrate the effectivity of our algorithm. 
Algorithm 7 The inputs are f,g G Q[x, y,z] such 
that gcd(/, g) = 1 and satisfying the two assump- 
tions, a bounding box B = x [Yi,!^] X \Z\,Z^\ 
and an error bound e > 0. The outputs are piece- 
wise rational parametric regular curve segments C, := 
{(x,y l (x),z i (x))( or (x l (y),y,z i (y)),a i < x(ory) < 
bi, (i = 1, N)}, which give a G 1 -continuous approxima- 
tion to f A g m B with precision e. 

(i) Topology determination and regular curve segmenta- 
tion of the plane curve defined by C\ : n z (f A g). 

(ii) Compute a rational number s as mentioned in Theo- 
rem 4. 

(iii) Let ip s : [x, y, z) — > (x,y + s z, z). Topology determi- 
nation and regular curve segmentation of the plane 
curve defined by C 2 : n z (ip s (f) A (p s (g)). 

(iv) Find out the correspondence between the regular 
curve segments of C\ and C2 ■ 

(v) Approximate the regular curve segments without VT 
point of C\ and C2 with eo < / f +4 e an d the ones 
with VT point with precision eo < 2 ^ +4 e - 

(vi) Recover the space approximation regular curve seg- 
ments of / A g with formula (10). 

(vii) Re-parameterize the non-rational approximation 
curves to rational approximation curves under the 
error control if there exist. 

(viii) Output the piecewise approximation regular curve 
segments. 

We will show several examples to illustrate our algorithm. 



Example 1 Consider the algebraic space curve defined by 
the system{f,g} = {x 2 +y 2 + z 2 -A, (z-1) {x 2 +y 2 -3 z 2 )} . 
In fact, they are two plane circles with z = ±1 as shown in 
Figure 7 (green ones). The space curve is not irreducible, 
not regular, and not in a generic position. We will approxi- 
mate it with rational curves under precision 1CP 2 . Following 
Algorithm 7, we have 

(i) Compute the resultant of f, g w.r.t. z, we haveC\ : h = 
x 2 +y 2 —3 = 0, as the red circle in Figure 7. We splitCi 
into eight regular curve segments with x- coordinates 
[-1.732050808,-1.0,0,1.0,1.732050808]: And the 
tangent directions of the points all are (1,— f^/f^) 
evaluated at these points. Note that x — corre- 
sponds x- critical points of C± and (±1.732050808,0) 
correspond to VT points, 
(ii) Since 7Tj,(/i) = x 2 — 3 ; we can obtain 

a = -2,ai = -1.732050808, ot 2 = 0, q 3 = 1.732050808, ct 4 = 2. 

And we can get r — 3.464101616. Following Theorem 
4, we have R = 1.0 when we choose g to compute R. 
We can select s = l<^ R = 1.732050808. 

(iii) Compute the resultant of ip(f) = f(x,y + z,z),ip(g) = 
g{x,y + z,z) w.r.t. z, we have C2 : h = (x 2 + y 2 — 
2 + 2 y) (—2 + x 2 — 2 y + y 2 ) = 0, as two blue circles 
in Figure 7. Since Tr y (h) — (x 2 — 3) (x 2 — 2) = 0, we 
split C2 into 16 regular curve segments at 

x = {-1.732050808,-1.414213562,0, 1.414213562,1.732050808}. 

And the tangent directions at the endpqints of 
these regular curve segments are (l>fj/§^) eval- 
uated at the points close to these points. We can 
get the approximating tangent directions. And we 
can find that (±1.732050808, ±1) are VT points 
since the absolute values of evaluated at 

(±1.73204, ±1.006118660), (±1.73204, ±.9938813400) 
are larger than 200. 

(iv) As shown in Figure 7, the critical points of C\ are 
P,Q. Choose a vertical line which intersect C\ at 
W,T. K(K = P,Q,W,T) are points on d and 
K\,K2 are corresponding points of K on C2. Con- 
sider W(0, 1.732050808), Wi(0, 2.732050808), W 2 (0, 0.732050808) 

for example. We can find that W\ , W2 are on the 
line x = in a neighborhood with radius 1.732050808 
centered at W. So we can conclude that W\, W% cor- 
respond to W with local generic position method. The 
correspondence of other points are similar. 

(v) Approximate C\,C2 respectively. In order to derive 
the required precision 10 _2 ; we we use precision e\ = 
0.0044 < Vl ^ +4 10~ 2 for the regular curve segments of 
C\,C2 without VT point(s), and we use precision e 2 



Fig. 7. Projection curves and approximation curve in Example 1. 



0.0022 < 2v / 1 1 2+4 10~ 2 for the regular curve segments 
with VT point. Consider a regular curve segments on 
Ci, (-1.732050808,0), (-1.60, 0.6633249580) are the 
endpoints for the one, denoted as C\ . And it has a VT 
point. (-1.60, 0.6633249580), (-1.40, 1.019803903) 
are endpoints for the other, denoted as G%. And 
it has no VT point. The approximation of C\ 

is {x, 0.9999999059 \J -x 2 + 0.000000464 x + 3.000000806), x S 

[-1.732050808, -i.6o] and the error is very small. The 
approximation for C 2 is (x, 0.6106757885 x + 2.310809554 - 

0.1270414345/(0.5070598449 x -f- 1.0), x e [-1.60,-1.40] and the 

error is 0.0004 < ei. For the regular curve seg- 
ments on C 2 with endpoints: [(—1.732050808, 1.0), 
(-1.60,1.663324958)], denoted as C 3 and it has a VT 
point. (-1.60, 1.663324958), (-1.414213562,2.0) are 
endpoints for the other, denoted as C4, without VT 
point. Similarly as C\, the approximation for C3 

is (x, 1.0 + 0.9999999054 */ -x 2 + 0.000000466 x + 3.000000809, x 6 

[-1.732050808,-1.60]. The approximation for C4 is 

(x, 0.6301674345 x + 3.324059999 - 0.1265242054/(0.5086364591 2! + 

1.0)), x g [-1.60, -1.414213562] and the error is 0.0002 < e\ . 
We can find that parts of C\ , C 2 and C3 , C4 are 
correspondent. 

(vi) Recover the approximation space curves of f A g by 
the formula z — V 2 ^^ 1 ^ , The space regular curve 
segment corresponding to C\ and C3, we have its ap- 
proximation parametric space regular curve segment 
for x £ [-1.732050808, -1.60]: 



(x, 0.9999999059 \j -x 2 + 0.000000464 x + 3.000000806, 
1.0 + 0.9999999054 -x 2 + 0.000000466 x + 3.000000809 
-0.9999999059 \J -x 2 + 0.000000464 x + 3.000000806). 

The approximation space curve is not rational, de- 
noted as Si . The approximation corresponding to 
C 2 ,C4forxe [-1.60,-1.414213562], denoted as S 2 , 
is 

(^,0.6106757885 x + 2.310809554 - .5 705 984 49 3 *+1.0 ■ 

0194916460 x + 1 013250115 0.1265242054 , 0.1270414345 \ 

U.U1J4J1D4DU Itl.lllMJU«J 0.5086364591 -c + 1.0 + 0.5070598449 i + l.o'' 

(vii) We will re-parameterize S\ into rational one. At 
first, we can find that the y coordinate of Si changes 
from to 0.6633249580. Its two endpoints are 
P (-1.732050808, 0, l.Q), P x (-1.60, 0.6633249580, 1.0) 
The tangent direction ofS 2 atP 2 is (1, 2.412090757, 0.0) 
By approximating the tangent direction of Si at Pi , 
we have (1,283.0783218,0.0). And there is another 
regular curve segment which shares a same tangent 
direction with Si at Pi . Taking their average value, we 
can set the tangent direction of Si at Pi as (0, 1,0). 
Using Formula (11), we can easily obtain the rational 
approximation regular curve segment for Si is 

(-2.412090984 y - 22.71853197 + _ llfgfs? li g£ y + 1 , V, 1-0). 

The error in x-direction is 0.0020563160 < e 2 , (We 
take 19 sample points besides endpoints to compute 
the error.). So the approximation rational curve sat- 
isfies the error requirement, 
(viii) Output the piecewise approximation curve. 



Example 2 Approximate the algebraic space curve de- 
fined by f — g = 0, where f = x 2 + y 2 + z 2 — 4, g = 
(a; 2 + y 2 + 2 y — z 2 ) (2 — x — 4y). It is a space curve with 
singular point. The approximation space curve is as the left 
part of Figure 8 and the error is 0.013. The color differs 
the different approximating space regular curve segments. 




Fig. 8. Appriximation curve and projection curves in Example 2. 

Example 3 In this example, we will approximate the al- 
gebraic space curve defined by f = g = inside [—2, 2] x 
[-2, 2] X [-2, 2] with errors = 0.014, where f = 95-94 x 3 + 
64 x 2 y + 28 x 2 z - 61 x 2 + 69 xy 2 - 53 xyz -59xy + 28 xz 2 - 
15 xz - 83 x - 3 y 3 + 59 y 2 z + 49 y 2 + 4 yz 2 + 1 1 yz + 5 y - 
81 z 3 - 8 z 2 - 9 z, g = 49 + 7 x 3 - 46 x 2 y + 87 x 2 z + 94 x 2 + 
73 xy 2 + 93 xyz - 3 xy - 27 xz 2 + 56 xz + 70 x + 72 y 3 - 
37 y 2 z -20y 2 + 79 yz 2 - 78 yz - 3 y + 94 z 3 + 30 z 2 + 47 z. 
The approximation space curve is as Figure 9. 




Fig. 9. Appriximation curve and projection curves in Example 3. 

5. Conclusion 

We introduce a local generic position method to com- 
pute the topology as well as the piecewise approximation 
curves of algebraic space curves. Especially, we present an 
algorithm to approximate algebraic space curve by piece- 
wise rational curves with correct topology and under any 
given precision. The method is effective. 
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